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Abstract 

A plasma becomes quantum when the quantum nature of its particles significantly affects its 
macroscopic properties. To answer the question of when the collective quantum plasma effects are 
important, a proper description of such effects is necessary. We consider here the most common 
methods of description of quantum plasma, along with the related assumptions and applicability 
limits. In particular, we analyze in detail the hydrodynamic description of quantum plasma, as 
well as discuss some kinetic features of analytic properties of linear dielectric response function in 
quantum plasma. We point out the most important, in our view, fundamental problems occurring 
already in the linear approximation and requiring further investigation. (Submitted to Physics- 
Uspekhi) 
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I. INTRODUCTION 



Plasma can be regarded as quantum when the quantum nature of its particles signifi- 
cantly affects its macroscopic properties. To determine when quantum effects in plasmas 
are important, their adequate description is needed. As plasma is an ensemble of many 
particles, the corresponding approach must be based on the appropriate description of a 
quantum particle. 

The description of non-relativistic quantum plasma can be based on either Schrodinger's 
representation (in which the operators are time- independent, while the time dependence of 
physical quantities of the system is defined by the corresponding time dependence of the 
system's wave function or density matrix) or Heisenberg's representation (in which the time 
dependence is transferred from the wave functions to the operators). Most of quantum 
plasma models that are commonly used now |l] use the Schrodinger's representation; the 
quantum plasma state is described either by wave functions of separate particles (the so- 
called multistream model {2!), or by the density matrix, or by the Wigner function defined 
in terms of the density matrix in coordinate representation , or - and this approach 
has become popular recently - by a set of the so-called quantum hydrodynamics equations. 
Naturally, simplifying assumptions are made in all these models (which thus lead to limi- 
tations of their applicability), which one should take into account when analyzing results 
obtained from them. However, concrete applicability limits of results obtained from a par- 
ticular model are not always stated explicitly (this especially concerns the widely used model 
of quantum hydrodynamics {5,6]), which can lead to their incorrect interpretation. This has 



recently been pointed out, for example, by Melrose and Mushtaq [3] as well as by Kuzelev 
and Rukhadze 8|. 

With the recent rapid increase in the number of publications on quantum plasmas, the 
lack of detailed analysis of the made assumptions and the associated limitations for the 
most common quantum plasma models becomes increasingly obvious, and therefore it is 
useful to provide such analysis. Beyond this, there is an important problem of macroscopic 
observability of quantum phenomena in plasmas; this problem is also connected with the 
question of when (and which) quantum phenomena are important in quantum plasmas. 
The answer to this question of course depends on the models and approximations used to 
describe the quantum plasma. In this paper, we provide a detailed analysis of the quantum 



3 



hydrodynamics model, and study the kinetic features of analytical properties of the linear 
dielectric response function in isotropic unmagnetized quantum plasma. In doing this, we 
highlight the most important, in our view, fundamental problems associated with the linear 
response of quantum plasma, which require further investigation. 



II. BASIC METHODS OF DESCRIPTION OF QUANTUM PLASMAS 

As in the case of classical plasmas, the most complete description of quantum plasma as 
a system of many interacting particles is a completely hopeless task. In case of quantum 
plasma this task is in a sense even more hopeless than in case of classical plasma, not only 
because it is impossible to solve the Schrodinger's equation for the A^-particle wavefunc- 
tion of the system, but also because of the lack of such wavefunction for a macroscopic 
system that interacts, however weakly, with its environment [91]. Yet the problem can be 
significantly simplified by assuming that the plasma is nearly ideal, i.e., that the two- and 
higher-order correlations between its particles can be neglected. If this is the case, then the 
plasma can be considered as a collection of quantum particles that interact only via their 
collective field. As mentioned above, the most commonly used now l| are the following 
models (which all in fact use the assumption of ideal plasma): (i) the quantum analog of the 
multistream model [2], (ii) the kinetic model based on the Wigner equation for the density of 
quasi-probability of particle distribution in coordinates and momenta 3|, y], and, finally, (iii) 
the quantum hydrodynamics model. All of these models, in one way or another, are based 
on the Schrodigner's equations for the wavefunctions of plasma particles, and therefore are 
non-relativistic; hence they can only be used for describing non-relativistic ideal plasmas 
and, strictly speaking, for describing plasma oscillations with small (non-relativistic) phase 
velocities uj/k <^ (here u and k are frequency and wavenumber of the oscillations, 

respectively, and c is the speed of light). We should note, however, that more general rela- 
tivistic models of "quantum plasmadynamics" appeared recently 10|; we expect that these 
models will be more widely used in the future, due to their logically consistent description 
of both quantum particles and quantized fields. However, here we will only consider the 
non-relativistic models, as they are the most widely used in recent literature on quantum 
plasmas, probably owing to their relative simplicity. 

Description of quantum plasma should be started with the models of Hartree and Hartree- 
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Fock, in which independent Schrodinger equations for plasma particles are coupled via 
the average self-consistent field (in Hartree-Fock's model the correction due to exchange 
interactions is also taken into account). The Hartree and Hartree-Fock approximations for 
quantum plasmas are analogous to the self-consistent field approximation in description of 
classical plasmas, and they form a basis [lay a foundation?] for kinetic and hydrodynamic 
models of quantum plasmas. Therefore it is important for further considerations to write 
down their main assumptions. 

Main assumptions of the Hartree and Hartree-Fock models: 

1. Plasma particles interact only through average classical (i.e., not quantized) collective 
fields. 

2. Plasma is ideal, F^ = f/mt/cF = e^n^^'^/ep ~ (TiUp/tFY ^ where ep = 
(^V2"^)(37^2 n)^/^ is the Fermi energy of electrons, Up = (47re^n/m)^/^ is electron 
plasma frequency, and e, m and n are charge, mass, and number density of electrons. 
We should note that for the electron gas in metals the condition Fg ^ 1 is in general 
not satisfied: in metals we have F^ ~ 1. 

3. Non-relativistic approximation is used; see the discussion after Eq. (!23l) . 



Kinetic models of Wigner-Poisson and Wigner-Maxwell These models are based on 



the Wigner equation describing time evolution of the Wigner function {4], which is coupled 
with either Poisson's equation or Maxwell's equations describing the self-consistent collective 
electrostatic and electromagnetic field, respectively. The Wigner function describes quasi- 
density of quantum particle probability distribution in coordinate-momentum phase space 
(we call it quasi-density because the Wigner function can attain negative values, due to 
noncommutativity of position and momentum operators in quantum mechanics, i.e., due to 
uncertainty principle). The Wigner function /(q, P,t) is defined from the density matrix 
p(q, q',t) of plasma in coordinate representation as follows [3|: 

1 r ..^.p / 1. ^ 1 



/(q, P) = ^ J dre-'-'-p (^q - -hr, CL+'hrj, (1) 

where q and P are canonically conjugated generalized coordinate and momentum, is the 
number of components of P (and/or of q) and corresponds to the number of coordinates of 



a particle (A^ = 3 in a 3-dimensional system). The Wigner function is normalized so that 

^(q)= [ /(q,P)rfP, (2) 



where n(q) is the number density of plasma particles. The Wigner equation, which describes 
evolution of the Wigner function, follows from the evolution equation for the density matrix 
in coordinate representation, and has the following form [sl. 



dt {27r)^h 



X 



(3) 



H { ff-\ — hk, r hr] — H in /ik, r H — hf 

'' 2' 2/ I' 2 ' 2 



where H is the Hamiltonian of the system. 

For a system of charged particles interacting via a self-consistant electrostatic field with 
a potential 0(q), the Hamiltonian is if(q, P) = p^/2m + e0(q) = i^(q, p), where p is the 
kinetic momentum of a particle (which in this case coincides with the generalized momentum 
P). For such systems, the Wigner equation ([3]) becomes 



(4) 



Equation (jl]) coupled with the Poisson's equation for the electrostatic potential 0(q) (in 
which the density of charged particles (electrons) is defined by Eq. ([2])), describes a weakly 
correlated system of charged particles interacting electrostatically, and is called the Wigner- 
Poisson model. 

For a system of spinless charged particles interacting via a self-consistent electromagnetic 
field, the Hamiltonian is i/(q, P) = (P — eA(q)/c)^ /2m + e0(q), where P is the canonical 
momentum of a particle, and 0(q) and A(q) are scalar and vector potentials of the electro- 
magnetic field, respectively. By changing variables according to P = p + eA(q)/c, where p 
is the particle kinetic momentum, the Wigner equation ([3]) is then cast in the form (we note 
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that the corresponding equation (30) of Ref. [3] has typos in some of the signs 
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(5) 



Equation ([5]), coupled with the Maxwell's equations for the self-consistent electromagnetic 
field and a gauge condition (e.g., the Coulomb gauge V-A = 0), describes a weakly correlated 
system of charged spinless particles interacting electromagnetically, and is called the Wigner- 
Maxwell model. 

Main assumptions of Wigner-Poisson and Wigner-Maxwell models: 

1. Plasma is ideal, Tq = t/mt/eF = e^n^/^/tp ~ {hL/jp/tpY ^ 1- As noted above, this 
condition is not satisfied for electron gas in metals, where Fg ~ 1. 

2. Plasma particles interact only via average collective fields that are described by 
Maxwell's equations (i.e., classical electrodynamics is assumed for the fields). 

3. Collisions between quantum particles are not taken into account (the models are col- 
lisionless) . 



4. Non-relativistic approximation is used; see the discussion below after fl23|) . 

5. Usually, spin of particles (as well as exchange interactions) are not taken into account. 
However, the effect of spin can still be accounted for in the workframe of non-relativistic 
Wigner-Maxwell model by introducing the spin distribution function and writing the 
corresponding kinetic equation for this function; this has been done, for example, by 



Silin and Rukhadze 



12| 



Multistream model 

This model is based on the Hartree approximation of plasma particles interacting via 
self-consistent collective fields only. Plasma is considered as a collection of weakly corre- 
lated "cold beams" formed by groups of particles with the same momenta; these beams 
are assumed to interact only through their collective fields. Using linearized equations of 
"cold hydrodynamics" for each of these groups of particles (beams), their current density 
and the corresponding dielectric permittivity tensor are calculated. Then adding up the 
contributions of all groups of plasma particles with the corresponding "weight functions", 
i.e., averaging these contributions over plasma equilibrium distribution function /o(p), one 
obtains the dielectric permittivity tensor of the whole plasma (see Sec. IIVI for details of this 
procedure). This procedure is equivalent to calculating the dielectric permittivity tensor 
directly from the Wigner equation ([5]), since the latter is based on the same assumption of 
weakly correlated particles interacting only via their collective fields. 

Main assumptions of the multistream model: Same as for kinetic models based on 
the Wigner equation (see above). 

Model of quantum hydrodynamics 

This model, first described in Ref. 13|, is constructed similar to the multistream model 



(note that it can also be derived from the Wigner model [14|]): the wavefunctions of plasma 



particles are represented in the form [l5| '?/'Q,(r, t) = aa{r,t) exp{iSa{'r,t)/h), where aQ,(r, t) 
and Sa{r,t) are the real functions of space and time, a is the particle index. The density 
Ua and velocity Vq, of an a-th particle are defined as Ua = |^/'Q,(r,t)p = a^(r, t), Vq, = 
VSa{r,t)/'m. The macroscopic plasma density n{r,t) = (jia) and velocity u(r,t) = (vq,) 
are introduced, where (...) denotes averaging over an ensemble of plasma particles, and two 
equations are written for n{r,t) and u(r,t): the equation of continuity and the equation 
of motion, the latter of which contains two pressure-like terms [l^ - the classical pressure 
defined as P'^^ = mn{{vD — (fa)^), and the quantum pressure 

^' = 7^((Va„)'-a„(VV)). 

In order to close the set of these two equations, the following two assumptions are made 
for P« and P^': 
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1. It is assumed that wavefunctions of all plasma electrons have equal amplitudes 
aa(r, t) = a(r, t) (which nevertheless can vary in space and time), while having different 
phases ^^(r,^). This assumption is in agreement with the assumption of uncorrelated 
plasma particles: indeed, the spatial distribution of each quantum particle, defined 
by the amplitude aa{r,t), is independent of the spatial distribution of other particles 
in the system. The assumption aa(r,t) = a{r,t) implies that the spatial distribution 
density of each quantum particle TZq/ — I GjQ^ I IS proportional to the density n of the 
whole system of particles, i.e., all particles are "smeared" over the whole system in 
the same way; in other words, the size of the wavepacket representing each particle of 
the system is equal to the size of the whole system of particles. (We note that this 
assumption does not impose any restrictions on the spatial and temporal scales of the 
waves that can be correctly described within this model.) This assumption implies 
the following relation between P'^ and n [14|]: 



2. Some equation of state is assumed, which links the "classical" pressure of the quantum 
plasma P^' = mn((f^) — (fa)^) with the macroscopic density of plasma n(r, t) = {ua)- 

Assuming a particular equation of state for the classical pressure P'^' imposes the cor- 
responding limitations on the applicability of thus obtained hydrodynamics model; this is 
discussed in more detail below in Sec. IIIII Beside the assumptions 1 and 2 above, the hy- 
drodynamics model of quantum plasma also makes all the main assumptions of the kinetic 
models based on the Wigner equation (see above). 

To summarize, we list all the main assumptions of the hydrodynamics model of 
quantum plasma: 

1. Plasma is ideal, Tg = U-mt/^F = e^n^^^/ep ~ {hiOp/epY ^ 1- As already noted, this 
condition is not satisfied for the gas of conduction electrons in metals, where Fg ~ 1. 

2. Plasma particles interact only via average collective fields that are described by 
Maxwell's equations. 

3. Collisions between quantum particles are not taken into account (the models are col- 





lisionless) . 
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4. Exchange interactions between plasma particles are ignored. 

5. Non-relativistic approximation is used; see the discussion below after (12^ . 

6. The wavefunctions of all plasma electrons ipa{T',t) = aa{r,t) exp{iSa{r,t)/h) are as- 
sumed to have equal amplitudes aa{r,t) = a{r,t) (which nevertheless can vary in 
space and time), while differing in their phases ^^(r,^). This imposes a relation be- 
tween P'^ and n - the "equation of state" ([6]) for the quantum pressure P'^. Note that 
this assumption restrains the size of plasma particle wavepackets (requiring them to 
be equal to the size of the whole system), but does not restrain the frequencies and 
wavenumbers of plasma waves that can be adequately described by this model. 

7. Some equation of state is assumed, which relates the "classical" plasma pressure P*^' = 
mn{{vD — (ffv)^) with the macroscopic density of plasma n(r, t) = {ua)- Usually the 
adiabatic equation of state is postulated, P^^ = PQ^n/rio)^, with Pq = riQep for 
degenerate electrons (i.e., for Tg ^ e^?, where Tg is the electron temperature in energy 
units), or with Pq = noTg for non-degenerate electrons (i.e., for Tg ^ ep). From 
this equation of state follows the following restriction for waves that can be correctly 
described within such hydrodynamics model: kXp 1 for degenerate electrons, where 
Xp = vp/VSujp is the Fermi- Thomas length,, vp = ^/2ep/m is the Fermi velocity of 
electrons, or /cAd ^ 1 for non-degenerate electrons, where Xd = A/Te/27re^n is the 
electron Debye length (see Sec. IIIII below). 

III. ON APPLICABILITY RANGE OF QUANTUM HYDRODYNAMICS EQUA- 
TIONS 

In some works (see the reviews by Manfredi et al. p. 26, after Eq. (54), and by 
Manfredi p. 14, the discussion after Eq. (4.30)) the following statement is made: "it can 
be shown that, for distances larger than the Thomas- Fermi screening length Lp {Xp in our 
notations), one can replace Ua with n" in the quantum pressure 
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As a result of this replacement, the set of hydrodynamics equations becomes (in one- 
dimensional case) 

du du _ e dcj) 1 dP"^ 1 dP'' 
dt dx m dx mn dx mn dx ' 

where P'^ = P'^{n) is defined by Eq. ([6]). This statement in fact implies that replacing 

with n = (ua) in (JTj), which is equivalent to postulating the "equation of state" ([6]) for 

the quantum pressure, is only valid at length scale large compared to Xp- This essentially 

means that postulating the "equation of state" ([6]) for the quantum pressure P'^ imposes 

the condition that the lengths of waves correctly described by Eqs ([8])-(|9]) should be large 



compared to the Thomas-Fermi length: kXp <C 1 



ep). The proof of this statement in Refs 16| and 



"or degenerate electrons, when ^ 
l|) is based merely on the fact that 



only for /cA^ ^ 1 the hydrodynamics equations (jH])-© correctly describe the dispersion 
of longitudinal oscillations in degenerate electron gas (this follows from comparison of the 



dispersion relations derived hydrodynamically and kinetically, which is done in Refs [16| and 
hence the limitation kXp 1 must be imposed somewhere during the derivation of 
the closed set of hydrodynamics equations (EJ-dH]) (which is correct), namely (and this is 
incorrect!) - when the "equation of state" (E]) for the quantum pressure P'^ is postulated, 
which is equivalent to replacing ria with n in ((Tj). 

It seems appropriate to clarify this issue here. Let us show that, for plasma with degener- 
ate electrons (i.e., for <^ ep), the limitation kXp 1 of quantum hydrodynamics appears 
not as a consequence of postulating the "equation of state" (|6]) for the quantum pressure P'^, 
but as a consequence of postulating a particular equation of state for the classical pressure 
pel _ pci^^'^^ namely - the adiabatic equation of state P^' = PQ''{n/noY , with Pq' = rioep 
(where uq is the equilibrium electron density); it is this equation of state that leads to the 
following dispersion of longitudinal oscillations in degenerate plasma 6| : 

.^ = .,^ + 5fc^4 + (i + «)|^ (10) 

where K = {'i8/175)mh^p/h^uj^. 

Without any additional assumptions, except the assumption of ideal nonrelativistic 
plasma, the equations for the plasma density n{x,t) and hydrodynamic velocity u{x,t) 
in one- dimensional case have the form ([8])-([9]), with the "quantum pressure" P'^ = P'^{na) 
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defined by Eq. ([Tj), and the "classical pressure" P'^^ defined as P'^'(r, t) = mn{{v^) — (fo)^) (it 
is called "classical" because it corresponds to the gas pressure in the classical limit h ^ 0, 
while the "quantum" pressure P'^, from which the Bohm diffusion term in appears, 
does not have an analog in classical plasma). As already noted, in order to close this set of 
equations, one needs to introduce two simplifying assumptions: (1) postulate an equation 
of state for the classical pressure - a relation between P'^' and n, and (2) postulate the 
"equation of state" (|6]) for the quantum pressure P"^, which is equivalent to replacing Ua 
with n in ([7]). Let us consider these two assumptions separately and see which limitations 
they impose on the applicability of the resulting set of hydrodynamics equations. 

First, we consider the assumptions implied by postulating the equation of state for the 
classical pressure; to do this, we consider the classical limit of Eqs OH])-® (to exclude, for 
now, the Bohm diffusion term with P'^, which vanishes in the classical limit). We note 
that in classical plasma nothing prevents us from considering the Fermi-Dirac distribution 
of electrons as the equilibrium distribution - in the classical plasma we are free to construct 
any equilibrium distribution at our will (as long as it is stable). As indicated in the textbook 
by Aleksandrov et al. in collisionless plasma (for which Eqs (I8])-(|9]) are written) there 
are two cases when the pressure can be evaluated directly, and the set of hydrodynamics 
equations ([8])-([9]) can be closed. The first case corresponds to processes with characteristic 
lengths L and times r, whose characteristic velocity greatly exceeds either electron thermal 
velocity (for Maxwellian distribution of electrons) or electron Fermi velocity (for Fermi-Dirac 
distribution of electrons): 



In this case, following 



— ~ — > max|fT,fF|- (11) 
r k 



17| , we can completely neglect the thermal or Fermi spread of electron 
velocities, which then corresponds to the case of cold plasma (ions are assumed to be cold), 
and from (j9]) obtain the Euler equation with zero pressure, P^' = 0. Naturally, this approx- 
imation does not yield the correction to the dispersion of plasma waves due to thermal or 
Fermi velocity spread of electrons; in order to obtain this correction, one needs to take into 
account thermal or Fermi velocity spread of electrons, which is done further below. 
The second case corresponds to processes for which 

L a; , , 

^^Ti < — ~ T ^ max|fT,t^F|- (12) 
T k 

In this case the effect of electron inertia is negligibly small (so that electrons have Boltzmann 
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distribution), and by excluding the electric field from the momentum equations for electrons 
and ions one obtains the set of one-fluid hydrodynamics equations. These equations are 
suitable for description of processes such as ion sound (with the hmitation (TT^). but are not 
suitable for description of electron oscillations, for which the electron inertia is essential; we 
will thus not consider this case here. 

Let us return to the first of the two cases mentioned above - the case of fast processes ffTTj) 
- and take into account the effect of velocity spread of plasma electrons (the ions are still 
regarded as cold). We consider small perturbations of equilibrium that is characterized by 
the Fermi-Dirac distribution of electrons (emulating degenerate electron gas), so that the 
pressure at equilibrium is Pq' = noep (dropping the factor of order of unity). We write the 
equation for electron energy density (the 2-nd moment of the electron distribution function 
/) in one-dimensional case as 

?f^u?f,SP'-p^2'^^0, (13) 

CM OX OX OX 

where Q is the energy density flux defined asQ = (m/2) J {v ~ uf fdv. With the assumption 
of fast processes ffTTj) . the term with the energy density flux dQ/dx is small compared to the 
term dP'^^/dt, and can be neglected (hence the condition (ITTj) implies that the corresponding 
process is adiabatic). As a result, Eq. f|T3|) together with the continuity equation ([8]) becomes 

dt dx J ' 

from which follows the equation of state of electron gas for adiabatically fast processes 
pd jj^i =constant, or P^^ = PQ{n/nQ)^ with P^ ~ n-Qep. Substituting P^'^ = PQ{n/nQY 
in (in the classical case, i.e., without the Bohm term), we obtain the momentum equation 
for electrons that accounts for electron velocity spread in equilibrium (unlike the approxi- 



17j). It is this equation that yields the correction ~ k'^Vp (or 



mation of cold electrons in 
~ k'^Vp for non-degenerate electrons) to the dispersion of longitudinal electron oscillations. 

Therefore, the momentum equation of quantum hydrodynamics iQ with P'^' = 
PQ^n/no)^, Pq ~ riQep, which yields the dispersion of longitudinal electron oscillations, 
is obtained from the kinetic theory in the approximation of adiabatically fast processes with 
uj ^ kvp- With application to electron oscillations with uj ~ Wp, this approximation is valid 
for long wavelengths: /cAj? ^ 1 for degenerate electrons, oi kXo <^ 1 for non-degenerate 
electrons. 
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Let us go back to the quantum case. In this case: (i) degeneracy is no longer a mere 
result of construction of electron distribution function, but is an effect of quantum statistics 
(due to Pauli's exclusion principle), and (ii) the following term (the Bohm term) appears in 
Eq. (Q: 

2m dx ^-^ \Pn^ 

a * 

(the remaining terms of Eq. ([9]) are the same as in the classical case, and for them the 
discussion of the previous paragraphs can be repeated). For the Bohm term the following 
assumption is made: 

J^d_y^ d^^/dx^ _ d r d^y^/dx^ \ 
2m dx ^ \fnZ 2m dx \ \fn I 

Q, V \ V / 

i.e., ria is simply replaced by n in the Bohm term (hence, the "equation of state" for the 
quantum pressure is postulated). It can be easily seen upon substituting ria = iV'aP with 
ipaixjt) = aa{x,t) exp{iSa{x,t)/h) in the left hand side of ( 1T5]) that the assumption (113]) is 
equivalent to the assumption [6] of quantum hydrodynamics (see the list of assumptions at 
the end of Sec. [Tll), which by itself does not impose any constraints on the lengths of waves 
that can be described within the framework of quantum hydrodynamics, as already noted 
in Sec. ITTl 

Therefore, the equations of quantum hydrodynamics (IHI)-(l9]) with P^' = PQ^n/no)^, 
Pq ~ ^Qep, and P'^{n) defined by Eq. ([6]), used in a number of works considering longi- 
tudinal electron oscillations in ideal quantum plasmas, can be obtained from the kinetic 
theory in approximation of adiabatically fast processes (ITTl) . which for plasma oscillations in 
degenerate electron plasma is equivalent to the condition k\p <^1. It is important to stress 
that this condition occurs not as a limitation of validity of ( IT5l) (i.e., of the "equation of 
state" ([6]) for quantum pressure), but as a consequence of the approximation of adiabatically 
fast processes, for which the equation of state for the classical electron gas pressure has the 
form P^' = Pg'(n/no)^ with Pq ~ fiQep (in one-dimension case). In case of longitudinal 
oscillations in degenerate electron plasma, the condition fcA^ <^ 1 corresponds to the Lang- 
muir part of the spectrum u = {cOp + 3k'^Vp/5Y^'^, which is correctly described by quantum 
hydrodynamics, unlike the essentially kinetic part of the spectrum Up/vp ^ A; ^ mvp/h 
that is associated with (kinetic) resonance u ~ kvp + hk"^ /2m and is analogous to zero sound 
in an almost ideal Fermi g Sec. |V]for more details. 



14 



IV. DIELECTRIC PERMITTIVITY TENSOR OF QUANTUM PLASMA 



Linear response of quantum plasma to electromagnetic disturbance can be described 
by dielectric permittivity tensor eij(a;,k) of the plasma. It can be calculated by solving 
equations governing plasma dynamics in the presence of self- consistent electromagnetic field, 
linearized assuming that the fields and the perturbations of the medium by these fields is 
small. For collisionless quantum plasma, such equation is the Wigner equation ([5]). The 
procedure of calculation of eij{u,k.) from Eq. is rather cumbersome, however, one can 
obtain the same result in a somewhat simpler way, based on the quantum multistream model 
mentioned above in Sec. [TTl 

It is well known (see, e.g., Ref. [l^) that the linear dielectric permittivity tensor of 
classical plasma can be calculated, rather simply, using the classical multistream model. In 
this model, plasma is considered as a collection of uncorrelated groups of classical particles 
with definite momenta (i.e., cold beams), and for each of these groups the current density and 
the corresponding dielectric permittivity tensor of the group is calculated. The contributions 
of all groups of particles are then multiplied by the momentum distribution function of 
plasma, /o(p), and summed over, thus yielding the dielectric permittivity tensor of the whole 
plasma. This procedure is equivalent to calculating ejj(a;,k) directly from the linearized 
Vlasov equation with self-consistent electromagnetic field. 

Obviously, this procedure can be generalized to quantum plasmas as well, if the lat- 
ter can be adequately described as a collection of weakly correlated quantum particles, 
i.e., if the coupling parameter Tg ~ (hU p/ep Y is small. This generalization is done. 



for example, by Kuzelev and Rukhadze 



191]. Starting from the set of "cold" quan- 



tum hydrodynamics equations for a group a of quantum particles with wavefunctions 
4'a{T',t) = aa{r,t)exp {iSa{r,t) /h) (without accounting for spin, i.e., for spinless particles) 

dUa 



dt 



+ V ■ (n„v„) = 0, 



, e / V X B\ „ f 1 

+ (v, ■ V) v„ = - E + + -^V — 



dt m \ c J Aw? V n, 



g 

ja = eUaVa, V„ = VS'„ A, (16) 

c 

and linearizing them, it is easy to calculate the conductivity tensor (Tj"(a;,k) of the a-th 
group of particles, by expressing their current density via the electric field as jaj(u;,k) = 
afj{uj,'k)Ej{(jj^'k) (for the corresponding Fourier components). Tensor of dielectric permit- 
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tivity of the a-th group of particles is then obtained as 

4ni 



Ane^noa ( ^ , w - k ■ v« 



- k • v„ 



where riQa is the unperturbed number density of group a particles, ujk = hk'^/2m, and 
K = k/|k| is the unit vector along k. Summing up the contributions of all groups of plasma 
particles with corresponding unperturbed densities riQa = / dpafoa{Pa), the linear dielectric 
permittivity tensor of quantum plasma is obtained as 

(u — k ■ \Y ( - kiVj + kjVi k'^VjV^ 



2 f ' 



47re^ 

eij{uj,k) = 6ij - / rfp/o(p)- 



— k • v) 



2 1 



a; — k ■ V (w — k • v) 



iuj — k ■ 



■ (17) 



Writing Eijiu, k) of isotropic plasma in the form [17 1 



Sijiu), k) = e {u, k)KiKj + e ^{u, k){5ij — 



from (fT7|) we obtain the following expressions for longitudinal and transverse e**" dielectric 



permittivities of isotropic quantum plasma 



19|: 



e{u,k) = 1 + ^ I dp-— J-—, 



hk^ „ 



dp 



w — k ■ V 



^[/o(p)], 



(18) 
(19) 



where the difference operator £)[/o(p)] is defined as D[fo{p)] = /o(p + ^k/2) — /o(]9 — /ik/2), 
and Vj_ is the component of the particle velocity perpendicular to the vector k. Eqs f lTSj) 
and (fT9|) are equivalent to the expressions obtained directly from the linearized Wigner 
equation (j5]) for isotropic quantum plasma jl2l |: 

47re2 



e\u,k) = 1 + 
£*'-(a;,k) = 1 + 



uhk'^ 
27re2 



dp ^'J D[Mp)], 
cu — k ■ V 



J dp[k X v]2 I /^(p) + 



1 ^[/0(P)] 

huj — k- V 



(20) 
(21) 



where /o(p) is the derivative with respect to particle energy e = /2m. Indeed, it is easy to 
show that in isotropic plasma Eqs (!20|) and (!2T|) coincide with Eqs f|T8|) and (fT9|) . respectively. 
Besides, same expressions (|T8l) and (fT9l) are obtained for an isotropic quantum plasma by 
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Klimontovich and Silin [3|, and also by Kuz'menkov and Maksimov [20|] from their Eqs (27) 
and (28) for e\u, k) and 6^'^{u, k) in the hmit of no exchange interactions. 

Eqs fll8p - (IT^ are obtained using non-relativistic model of quantum plasma, which, strictly 
speaking, is not applicable for description of plasma waves with relativistic phase speeds, 
u/k ^ c, regardless of whether the plasma particle velocities are relativistic or not [7]. To 
see this, consider the processes of emission or absorption of a quantum (of longitudinal or 
transverse wave) with frequency u and wave vector k by a plasma particle with momentum 



p and energy e{p). In relativistic treatment, i.e., for e = ^Jvrfic^ + p^(? (where m is the rest 
mass of the particle, e.g., an electron), conservation of total energy and momentum in these 
processes 

e' = e ± /iw, p' = p ± /ik (22) 

leads to the resonance condition 

where 7 = (1 — v"^ / c^)~^/'^ . In case of nonrelativistic particle velocities, taking e ~ mc^(l + 
v^jlc?) p ^ mv(l + v^jlc?) in ([23]), we obtain 

However, if in nonrelativistic approximation one formally takes e = /2m in (!22|) . then the 
term uj^ j(? does not appear at all in the resonance condition ( 12^ . which, strictly speaking, 
is incorrect in case of relativistic phase velocity of a wave in resonance with the particle, i.e., 
for bjjk^c. And though the effect of the w^/c^ term turns out to be insignificant for most 
processes in unmagnetized plasma, it can be of crucial importance for some processes in the 
presence of an external magnetic field, e.g., for the process of cyclotron maser radiation |21 |. 
One should remember that, strictly speaking, the nonrelativistic approximation is justified 
a priori only when the phase speeds of waves (as well as velocities of plasma particles) are 
nonrelativistic [18|. Indeed, energy exchange in wave-particle interaction, Eq. ( l22i) . includes 
the possibility that, for example, an initially nonrelativistic particle may become relativistic 
after interacting with a wave quantum with sufficiently large energy/momentum. 

In relativistic treatment, the dielectric permittivity tensor of an isotropic electron- 
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positron gas (in which electrons and positrons are assumed unpolarized) has a form 

(25) 

where p = m7v, /o(p) = 2?2(p)/(27r^)^, n{p) is the sum of occupation numbers for electrons 
and positrons, and 



We note that the denominators of the integrands in the electron and positron contributions 
in (125|) correspond to products of the relativistic resonance conditions (!23|) for emission 
and absorption of wave quanta by relativistic particles (electrons or positrons). Indeed, 
symmetrizing ( 125|) on p under the condition of isotropic plasma, /o(— p) = /o(p)) and 
equating thus obtained denominator to zero 

[(c - k . v)2 - A^] [(a; + k.v)2-A^] =0, 

we obtain four conditions 

w±k-v = ±Afc, (27) 

which correspond to resonances in the processes of emission and absorption of photons 
or plasmons by electrons and positrons, as well as in the processes of one-photon or one- 
plasmon electron-positron pair creation and annihilation 10|. The longitudinal e\.^i{uj, k) and 



the transverse e^^^iiu dielectric permittivities of an isotropic plasma, obtained from (l25ll . 



match with the corresponding expressions obtained by Tsytovich 22 



typo in the signs in Ref. 



as adjusted for the 



10|, and also neglect- 



2^, corrected in Eqs (9.1.12)-(9.1.13) of Ref . 
ing the polarization effects). 

For electron plasma (without positrons) in nonrelativistic limit c — ?■ oo (i.e., for nonrela- 
tivistic plasma with 7 — 1, and for waves with u/k <^ c), Eq. (l25l) changes into 

^ir(^' k) = 5ij ^ / rfp/o(p)- 2 \ k + — 1 + 



(28) 



differing from 6ij{u, k) of f|T7|) in the term 



2 \ ^P/o P 7 T-\2 2 O^^S- - • 29 

muj^ ] (w — k ■ v)^ — 

The disagreement of the nonrelativistic limit ( |28|) of Eq. ( |25|) with the expression (fT71) ob- 
tained in several different ways from nonrelativistic models of quantum plasma is due to the 
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fact that Eq. ( 125|) is obtained for a gas of unpolarized particles (electrons and/or positrons) 
with spins ±1/2, while Eq. f lT7|) is obtained from the models that do not account for spin. 
We should note that the term f l29p . by which Eqs (128 p and (11 7p differ, contributes only 
to the transverse dielectric permittivity of plasma e**", and has no effect on the longitu- 
dinal dielectric permittivity eK Therefore, nonrelativistic theories not accounting for the 
spin are not quite correct: while correctly describing longitudinal plasma oscillations, they 
incorrectly describe transverse plasma modes. Note that accounting for pararnagnetic ef- 
fects (associated with the spin of plasma electrons) in the nonrelativistic model 12| leads to 
disappearing of the term f l2^ in f[T7|) : as a result, thus modified Eq. flTTl) becomes exactly 
equal to the nonrelativistic limit (125]) of the tensor (I25p . The discussed discrepancy between 
the nonrelativistic responses f[T7|) and fl25]) . associated with the spin of plasma particles, is 
another important example showing the necessity of careful consideration of all the relevant 
effects, however small they may seem a priori in the nonrelativistic approximation. 

To obtain the relativistic generalization of the collective linear response (fT7|) of a gas of 
spinless charged particles (e.g., a gas of Cooper electron pairs with zero total spins, or a gas 
of bosons with spin 0), one needs to construct the corresponding relativistic theory based 
on the Klein-Gordon equation (which is a special case of the Dirac equation for spinless 
particles, i.e., for spinors of the 1-st rank - scalars). However, one could instead employ 
a phenomenological approach to establishing the form of relativistic generalization of the 
dielectric tensor ( ITTll . Note that the only difference between the wanted relativistic response 
of a gas of spinless charged particles (which should turn into ( ITTl) when c -> oo) and the 
relativistic response of a gas of unpolarized electrons is in the term that turns into (l29l) in 
the limit c — )■ oo. The denominator of (l29l) in the nonrelativistic limit c — )■ oo corresponds 
to the product of two resonances fl23l) for emission and absorption of a wave quantum with 
energy hu and momentum hk. Therefore, in relativistic case this denominator should turn 
into the product of the corresponding relativistic resonances ( l23|) . i.e., ul should turn into 
A^, and (w — k ■ v)^ — should turn into (cu — k ■ v)^ — A^. Hence, in relativistic theory 
Eq. ( l29l) should turn into 

ine^ r A? 

muj^ J (w — k-v)^ — A^ 

Adding the term ( l30l) to ( |25|) . we obtain the following phenomenological expression for the 
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tensor of dielectric permittivity of relativistic quantum plasma of spinless charged particles: 



k 

c.-k-v + (a;-k-v)2 + (a;-k-v)2 /• ^"^^^ 
It is easy to verify that ( 13T|) turns into ( |T7I) in the limit c — )■ oo, as required. 

As discussed above, plasma response obtained from relativistic treatment (be it the re- 
sponse of unpolarized electron or electron-positron gas, or the response of a gas of spinless 
charged particles) differs from the corresponding plasma response obtained from nonrela- 
tivistic models, among other things, by a term proportional to u'^/c^. The effect of this term 
on the dispersion of longitudinal and transverse waves in nonrelativistic quantum plasma 
(with 7 = 1) is usually small in most cases (e.g., for electron gas in metals), but becomes 
significant at large plasma densities, when the process of pair creation by photons and/or 
plasmons becomes energetically allowed in plasma. In this case an additional mechanism of 
wave damping, associated with electron-positron pair creation, is switched on (in addition 
to Landau damping). This damping occurs for waves (both longitudinal and transverse) 



with superluminous phase velocities, and has the energy threshold 22 1 



{huY > A{mc'y + {nkYc\ (32) 

It follows that, in the limit of long wavelengths (when u ~ Wp), damping of longitudinal 
and transverse waves due to creation of real electron-positron pairs becomes significant at 



plasma densities n ^ 10^^ cm ^ 22|. This is, however, a rather exotic such large 

densities can exist perhaps only in the core of dense astrophysical objects, e.g., white dwarf 
stars. On the other hand, the effects associated with creation of virtual electron-positron 
pairs do not have the strict energy threshold (!32|) . and thus can affect the waves in plasmas 
with relatively low densities. They can also affect analytical properties of plasma linear 
response function. 



V. QUANTUM KINETIC EFFECTS AND ANALYTIC PROPERTIES OF LIN- 
EAR LONGITUDINAL PLASMA RESPONSE 

Let us now consider the essentially quantum effects that occur in the kinetic descrip- 
tion of collective modes in quantum plasmas. For simplicity, we only consider longitudinal 
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oscillations resulting from an initial perturbation of quantum plasma. 

In case of initial perturbation /(r, p, 0) (here / is the Wigner function, see Sec. HTj) of a 
uniform isotropic plasma, time evolution of the electrostatic potential in plasma is given by 
the following expression 

0(t,r) = — / 0^(u;,k)e-^"*du;, (33) 

where the integration is carried out in the plane of complex uj along the horizontal contour 
located in the upper half plane, Im(a;) = a > 0, and 

where is the component of particle momentum along k, (^(k, p^) = J (^(k, p)dpydpz, g(k, p) 
is the Fourier transform of the initial perturbation /(r, p, 0), and ^'(w, k) is the longitudinal 
dielectric permittivity of plasma, defined by Eq. (fTSll . 

To obtain the behavior of the potential (p at large times t, one needs to integrate over u in 
(!33|l along the contour formed from the initial contour by taking the limit a = lm.{u)) — > — oo, 
while preserving the analyticity of the function (pcui^^, k) under the integral. (In order to do 
so, the function </)(^(a;,k) in its turn should be analytically continued from the region of its 
definition Im(w) > to the region Im(a;) < 0; to do this, the contours of integration over p^ 
in the numerator and denominator of flM|) must be displaced from the real axis Im(p^) = 
into the lower half-plane Im(pj.) < in such a way that the pole p^ = mu/k is passed from 



below 23|].) This contour of integration over u in fl33|) must pass above all the singularities 
of the function (pcui^^, k) (analytically continued to the region lm{u) < 0) that lie on or under 
the real axis Re(ci;) |23j. In case of classical plasma, the equilibrium distribution function 
/o is an entire function of px (i.e., /o has no singularities at finite px), and the analytic 
continuation of £:'(u;,k) to the region lm{u)) < is also an entire function of u. The same 
argument applies to the analytic continuation of the function 

/ ; ^dpx, (35) 

in the numerator of Eq. (El]), if the initial perturbation 5f(k,p^.) is also an entire function of 
Px- Thus, provided fo{px) and g(k,px) are entire functions of px, the function 0(^(a;,k) (IMll 
is the ratio of two entire functions of u. Then the only singularities of (f)uj{^, k) are the poles 
defined by the roots of the equation e\u, k) = 0. The contribution of these poles into the 
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integral (133|) completely defines the evolution of r), which in this case is a superposition 
of oscillations, exponentially damping (or growing, in case of non-equilibrium plasma) with 
time. 

However, solution of the initial value problem in quantum plasma turns out to be more 
complicated ^], since the equilibrium distribution function /o is no longer an entire function 
of px- Indeed, an electron gas obeys Fermi statistics, and its equilibrium distribution function 
(Wigner function) is 

V/2m-/i(T)" 



^'^^^ = (2^ < 



T 



+ 1 , (36) 



where T is the temperature and /i(T) is the chemical potential of plasma electrons. In one- 
dimensional case {p = Px) this function (shown in Fig. [T]by dashed lines) has singular points 
(poles of first order) that lie in the complex p^ plane on hyperbolas intersecting the real 
axis \m.{px) = at points ±[2/i(T)/m]^/^ (see Fig. |2]), with the distance between adjacent 
singular points proportional to temperature T 2J], so that in the limit T — )■ these points 
completely fill the hyperbolas. As discussed above, when performing analytic continuation of 
the function e\uj, k) into the region Im(a;) < 0, the contour of integration over p^ in e\uj, k) 
must pass the pole Px = muj/k from below; however, this contour should not intersect any 
of the poles of fo{Px)- Therefore the contour of integration over px is "pinched" between 
the pole Px = mu/k (lying in the lower half-plane of complex px) and the nearest pole 
Pxo of the function fo{Px)- For mu/k — )■ pxo, i.e., when these two poles coincide, the 
contour unavoidably passes through (intersects) these two poles, and the function e^Ujli) 
thus has singularities at Uj = kpxoj/m, where kpxoj are the poles of fo{Px)- Thus in quantum 
plasma the function £:'(a;,k) is no longer an entire function of u, and the singularities of 
e''{u,'k) contribute into 0(t,r), along with the zeros of e\ijj,li). This contribution, unlike 
the contribution of the zeros of e^u, k), might be not of the form of an exponentially damped 
oscillating function, but rather can be a relatively slowly (by power law) decaying function 



of time 



24j . not necessarily oscillating. 



However, in most cases plasma electrons are described by three-dimensional, rather than 
one-dimensional, distribution function, fo{p) oc {exp[(p^/2m— /i(T))/T] + l}~^, where pis an 
absolute value of three-dimensional momentum. In this case the one- dimensional distribution 
function under the integral in (fTSl) is defined as the three-dimensional distribution function 
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Mpx) Mpx) 




three-dimensional (p = y Px + P^, shown with solid lines) Fermi-Dirac distribution (j36p , in case 
of complete degeneracy {^{T)/T — )• oo, left panel) and partial degeneracy {fi(T)/T = 10, right 
panel). The functions are normalized on their corresponding values /g^(0). 



Im(p,) 




FIG. 2: (Color online) Left panel: absolute value of one-dimensional Fermi-Dirac distribution 
fo{Px) {exp[(p^./2m — ^{T))/T] + 1}^^ as a function of complex p^. Right panel: poles of one- 
dimensional Fermi-Dirac distribution fo{Px) oc {exp[(p^/2m — ^{T))/T] + and the contour of 
integration over px (bold dashed line) for analytic continuation of e'(a;, k) into the region Im(a;) < 0. 
The hyperbole - the locus of the poles of fo{px) ~ is shown by dotted line. 

fo{p) integrated over momenta perpendicular to k: 

(it is shown in Fig. [T]by solid lines), while the difference operator Z)[/o(p)] = /o(p + /ik/2) — 
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/o(p - Kk/2) acts only on /o^^fe), ^[/o(p)] = [h{p. + nk/2) - fo{p. - hk/2)]fo{pi.). As a 
result, Eq. ( 1T8|) for e\ijj,k) becomes: 



e\u,k) = 1 + 



47re^ 47rmT 



dpx- 



In 



In 



1 + exp 



[u - kpx/m) 
{px - hk/2y - 2m/i(T)\ 



1 + exp 



2mT 



2mT 



(37) 



where the contour of integration over is again chosen to pass the pole Px = muj/k from 
below 231]. The logarithmic functions in (!37j) are not entire functions - each of them has 
singularities like branching points and branch cuts in the complex p^ plane, shown in Fig [3l 
These singularities lead to £:'(u;,k) having analogous singularities in the complex u plane, 
shown in Fig. HI which can contribute to the integral over uj in ( 133|) , along with the zeros of 
£'(a;,k). This contribution can also lead to a power-law, i.e., non-exponential character of 
temporal damping of an initial perturbation, similar to the case of one-dimensional Fermi 
distribution I24I . Note that this is a purely kinetic effect absent in quantum hydrodynamics 
model. 

Im(p,) Im(pJ 



branching 
points and cuts 
of logarittims 
in£'(a», k) 




Re(p>) 



brandling 
points and cuts 
of logarithms 
ine'(6;, k) 




Re(p,) 



FIG. 3: Singularities (branching points and cuts, shown with sohd hues) of the logarithmic functions 
under integral in (j37p . and contours of integration over px (shown with bold dashed lines) for 
analytic continuation of e'(a;,k) into the region Im(a;) < 0, in two cases illustrating possible 
location of the pole px = moj/k (shown with a dot) relative to the branch cuts. The dotted line 
shows the locus of the branching points (the hyperbole). 

If, under some conditions, the contribution of singularities of e\uj, k) of quantum plasma 
dominates over the contribution of zeros of e\uj,'k) into r) at large times t, then the 
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Im(a») 




Re(ai) 



FIG. 4: Singularities (branching points and cuts, shown with sohd hues) of e'(a;,k) in complex 
uj plane, and contour of integration over oj in (|33p . displaced into the lower half-plane to infinity 
and bypassing the singularities of e'(f^, k) that lie in the lower half-plane (shown with dashed line). 
Both the poles of the function l/e'((^, k) and its branch cuts in the lower half-plane of complex cj 
contribute to the integral in ([33ll . 



physical picture of longitudinal collective oscillation modes in quantum plasma is altered: 
in this case the evolution of an initial perturbation of quantum plasma at large times t 
can be qualitatively different compared to that of classical plasma. The question of when 
this can happen is of fundamental interest, as it is related to macroscopic observability 
of quantum effects in plasma, and requires further studies. Results of this study will be 
published elsewhere. 

We note, however, that even in cases when the contribution from integrating along the 
cuts in Fig. HJis small or zero, quantum kinetic effects still significantly affect the dispersion 
and damping of longitudinal oscillations. In case of completely degenerate (Te = 0) Fermi 
distribution of electrons, equation £:'(u;,k) = yields the following approximate solutions 
at long wavelengths (note a typo in the signs in the dispersion relation (1391) of the "kinetic 



mode" in Ref. 



25|): 



3 
5 



1/2 



^±{k) = (/cfi? ± cufc) {l + exp [— 2 — 2A;^A|^] } , Up/vp ^ k <^ mvp /Ti. 



(38) 

(39) 
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These limiting solutions are shown in Fig. |5l For k <^ Up/vp, the dispersion of longitudinal 
mode fl38|) corresponds to that of the usual Langmuir mode, which can also be obtained from 
the hydrodynamics model, see Eq. (ITU]) of Sec. Illlj thus we can call f l38l) the "hydrodynamic" 
mode. Yet for Up/vp ^ A; ^ mvp/h the dispersion of longitudinal mode f l5U]) is mainly 
defined by the kinetic resonances u — kvp ^ ooj. = Q (see Eq. (jM])) between plasmons and 
electrons whose velocities along k are equal to vp- For those values of k for which the 
"hydrodynamic" mode frequency (|38|) becomes close to the kinetic resonance frequency uj = 
kvF±u)k (which occurs near the intersection of the dispersion curve (138|) with the resonance 
line u = kvp + ook), kinetic effects become dominant and "switching" from the hydrodynamic 
mode ( 155]) to the kinetic mode ( 15U]) occurs. We also note that in degenerate plasma quantum 
recoil leads to an imaginary part of the longitudinal dielectric permittivity 2J] for k > ki ^ 
{uJp/vf)^/ (3/2) (I In ?7 1 — 1), where rj = hoOp/Aep ~ a/T^? which in turn leads to Landau 
damping of longitudinal oscillations with k > ki ^ {up/vp)^/{3J2)(J]nr]\'^^T) 24^] (while 
for k < ki Landau damping is zero). Thus quantum kinetic effects play an important 
(and for certain wavenumbers - dominant) role for dispersion and damping of longitudinal 
collective oscillations in quantum plasmas. 



























■ 

if 







FIG. 5: Dispersion of longitudinal oscillations of degenerate (Tg = 0) electron gas. The approximate 
solutions (|38p -p9 p are shown with dashed lines. The part of the dispersion curve for which Landau 
damping occurs {k> k\) is shown with dotted line. 
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VI. SUMMARY 



In this paper we analyzed applicability limitations following from the basic assumptions 
of quantum plasma models most widely used in recent literature. Lack of understanding of 
these limitations can lead to incorrect interpretation of results obtained with a particular 
model. For example, description of longitudinal oscillations of degenerate electron gas with 
the model of quantum hydrodynamics yields a dispersion relation that is only valid for 
kXp ^ 1, which is not always stated explicitly in the literature. Moreover, the very reason 
for this limitation on lengths of waves described by quantum hydrodynamics is not clearly 
understood by all. In Sec. IIIII we showed that the limitation kXp 1 appears in quantum 
hydrodynamics as a result of postulating a particular (adiabatic) equation of state for the 
"classical" pressure, and not as a result of postulating for the "quantum" pressure. 

We also discussed the linear response of quantum plasma, and pointed out that it is 
conceptually incorrect to use the response obtained from nonrelativistic models of plasma 
when describing waves with relativistic phase speeds u/k ^ c, regardless of whether the 
plasma itself is relativistic or not. Besides, we pointed out that, even for unmagnetized 
plasmas, ignoring the effect of spin leads to incorrect dielectric permittivity tensor, which 
in particular leads to incorrect dispersion relation for transverse waves. 

Finally, we discussed quantum kinetic effects associated with nontrivial analytic proper- 
ties of the complex linear plasma response function, occurring both due to quantum degener- 
acy of electron distribution and due to quantum recoil. These effects are hardly discussed in 
the literature, with rare exception, while being of fundamental interest. In particular, correct 
account of the analytic properties of quantum plasma linear response function can signifi- 
cantly change the physical picture of evolution of collective oscillations in both unbounded 
and bounded quantum plasma (Landau's initial value and boundary value problems 23|). 
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